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Abstract 

A tight correlation between 7-ray and radio emission is found for a sample of BL 
Lacertae (BL Lac) objects detected by Fermi Gamma-ray Space Telescope (Fermi) 
and the Energetic Gamma- Ray Experiment Telescope (EGRET). The 7-ray emission 
of BL Lac objects exhibits strong variability, and the detection rate of 7-ray BL Lac 
objects is low, which may be related to the 7-ray duty cycle of BL Lac objects. We 
estimate the 7-ray duty cycle, 5^ ~ 0.11, for BL Lac objects detected by EGRET and 
Fermi. Using the empirical relation of 7-ray emission with radio emission and the 
estimated 7-ray duty cycle 6-y, we derive the 7-ray luminosity function (LF) of BL 
Lac objects from their radio LF. Our derived 7-ray LF of BL Lac objects can almost 
reproduce that calculated with the recently released Fermi bright active galactic nuclei 
(AGN) sample. We find that ~ 45% of the extragalactic diffuse 7-ray background 
(EGRB) is contributed by BL Lac objects. Combining the estimate of the quasar 
contribution to the EGRB in the previous work, we find that ~ 77% of the EGRB is 
contributed by BL Lac objects and radio quasars. 

Key words: galaxies: active — galaxies: BL Lacertae objects: general — cos- 
mology: diffuse radiation — gamma rays: theory 

1. Introduction 

There are 66 high- confidence identifications of blazars in the third catalog of active 
galactic nuclei (AGN) detected by the Energetic Gamma-Ray Experiment Telescope (EGRET) 
on the Compton Gamma-Ray Observatory (GGRO), which include 20 BL Lac objects (e.g. 
Hartman et al. 1999; Mattox et al. 2001). The LAT bright AGN Sample (LBAS) detected 
by Fermi includes 104 blazars consisting of 42 BL Lac objects, 57 fiat spectrum radio quasars 
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(FSRQs), and 5 blazars with uncertain classification (Abdo et al. 2009a). Comptonization is 
believed to be responsible for the 7-ray emission from blazars, which can be classified into two 
categories: the external Comptonization (EC) model and the synchrotron self-Comptonization 
(SSC) model, according to the origin of the soft seed photons (see, e.g., Bottcher 2007 for 
a review and references therein ). In the EC models, the soft seed photons are assumed to 
originate from the accretion disks, the broad-line regions (BLRs), or/and the dust tori (e.g., 
Ghisellini & Madau 1996; Georganopoulos et al. 2001; Dermer & Schlickeiser 2002). It was 
suggested that the SSC model may be responsible for 7-ray radiation from BL Lac objects 
while the EC model predominates over SSC model for quasars (e.g., Dondi & Ghisellini 1995). 
The statistical analysis on a sample of EGRET blazars implied that the soft seed photons 
may be predominantly from the BLRs (e.g.. Fan & Cao 2004; Fan, Cao & Gu 2006), which 
is consistent with the lack of BLR and other external fields emission in most BL Lac objects. 
Their results are roughly consistent with the detailed modeling of SEDs for a large sample of 
7-ray bright AGNs by (Ghisellini et al. 2010), in which all external field photons from BLR, 
accretion disk, and dust torus, are properly considered. However, the situation becomes more 
complicated for radiation in other wavebands (e.g., see Kataoka & Stawarz 2005; Bai et al. 
2009; Liu et al. 2008). It was found that the X-ray radiation from some FSRQs (e.g., 3C 273) 
is consistent with combined SSC and EC mechanisms (e.g., Pacciani et al. 2009). 

The extragalactic diffuse 7-ray background (EGRB) was first discovered by the SAS 2 
satellite (Fichtel et al. 1978; Thompson & Fichtel 1982) and subsequently confirmed after the 
launch of EGRET (Michelson 1995; Sreekumar et al. 1998). The EGRB integrated above 100 
MeV was determined to be (1.45 ±0.05) x 10"^ photons cm~^ s~^ sr~^ (Sreekumar et al. 1998). 
Using a new model of the Galactic background. Strong et al. (2004) obtained a slightly smaller 
value of the EGRB, (1.14 ±0.12) x 10~^photons cm~^ s~^ sr~^. Almost all extragalactic 7-ray 
sources in the third EGRET catalog are identified as blazars, which account for about 13% of 
the EGRB (Hartman et al. 1999). The contribution of unresolved blazars to the EGRB has 
been explored in many previous works, either by extrapolating the observed 7-ray luminosity 
distribution to obtain a 7-ray luminosity function (LF) or using the correlation of 7-ray emission 
with the emission in radio bands derived from EGRET blazars, which showed that about ~25% 
to ~100% of the EGRB can be attributed to the unresolved blazars (e.g., Padovani et al. 
1993; Chiang et al. 1995; Stecker & Salamon 1996; Chiang & Mukherjee 1998; Miicke & Pohl 
2000; Narumoto & Totani 2006; Cao & Bai 2008). 

A new revised catalogue of EGRET 7-ray sources (EGR2008) was given by Casandjian 
& Grenier (2008), and a larger catalogue of the LAT Bright AGN Sample (LBAS) detected by 
Fermi was released recently (Abdo et al. 2009a). In this paper, we derive a correlation between 
7-ray and radio emission for a sample of BL Lac objects detected by EGRET or/and Fermi 
(Casandjian & Grenier 2008; Abdo et al. 2009a). We estimate the 7-ray duty cycle of BL 
Lac objects detected by EGRET and Fermi, and the 7-ray LF of BL Lac objects is calculated 
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Fig. 1. The relation of iT-corrected radio luminosity and 7-ray luminosity for a sample of 7-ray BL Lac 
objects detected by Fermi/EGRET. The solid line represents the linear regression for all the sources in 
the sample. Open circles: Fermi BL Lac objects; filled circles: EGRET BL Lac objects. 

with this empirical correlation from their radio LF and 7-ray duty cycle. The radio LF of BL 
Lac objects derived by Padovani et al. (2007) based on the Deep X-ray Radio Blazar Survey 
(DXRBS) is used in this work. We further estimate the total number of BL Lac objects to be 
detected by Fermi, and the contribution of all BL Lac objects to the EGRB. The cosmological 
parameters Qm = 0.3, = 0.7, and Hq = 70 km s~^ Mpc~^ have been adopted in this paper. 

2. The correlation between 7-ray and radio emission 

The recently released LAT Bright AGN Sample (LBAS) detected by Fermi includes 42 
BL Lac objects, of which 11 sources were already detected by EGRET. Till now, there are 51 
BL Lac objects detected either by EGRET or Fermi [with the Test Statistic (TS) value above 
the EGRET threshold (4cr for \b\> 10°, and 5cr for | |< 10°); or lOcr for the sources detected by 
Fermi]. There are 40 BL Lac objects with measured redshifts, which are adopted for our present 
investigation. The radio/7-ray flux and photon index F above lOOMeV for the sources are listed 
in Table 1. In figure 1, we plot the relation between radio luminosity Lr^sg at 5 GHz and 7-ray 
luminosity (> 100 MeV) for these BL Lac objects, where the average K-corrected radio and 
7-ray fluxes are adopted for calculating and Lr,5g- We find a significant correlation between 
Lr 5g and (the correlation coefficient r = 0.926). The linear regression of the sample gives 



Abdo et al. (2009a) investigated the relation between 8.4GHz radio flux density and 7-ray flux 
density for Fermi blazar sample. The peak values are adopted by Abdo et al. (2009a), which 
are likely to correspond to a short time-scale flare state of a source. Kovalev et al. (2009) 
investigated the relation between quasi-simultaneous 7-ray and 15GHz radio flux densities in 
Fermi blazar sample, which includes xx BL Lac objects. The mean fluxes above 100 MeV 



log i^L^ = 0.995 log i^Lr,5g + 2.886. 
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Fig. 2. The 7-ray LFs for BL Lac objects. The solid line is the LF derived in this work, while the dashed 
line and the dotted line are the LFs (redshift less than 0.3 and greater than 0.3) calculated with the LAT 
bright AGN sample detected by Fermi (Abdo et al. 2009a). 

are adopted in this work, and the sample is hmited to BL lac objects, which should be a 
suitable choice, because we focus on the statistic properties of 7-ray BL Lac objects as a whole 
population in this work. The significance of the correlation found in this work is higher than 
that found by Abdo et al. (2009a), which may be attributed to the fact that only the peak 7-ray 
fluxes were adopted in their analysis. For 7-ray quasars, the situation is more complicated, as 
their 7-ray emission is not only dependent on their radio emission, but also on the origins of 
the soft seed photons, which is beyond the scope of this work. 

3. The 7-ray luminosity function of BL Lac objects 

Based on the correlation between 50 and established with a sample of BL Lac 
objects, we can calculate the 7-ray LF of BL Lac objects from their radio LF. Recently, Padovani 
et al. (2007) derived a LF for BL Lac objects based on the Deep X-ray Radio Blazar Survey 
(DXRBS), which extends to lower luminosity than that derived by Stickel et al. (1991). In this 
work, we use Padovani et al. (2007) 's radio LF to derive the 7-ray LF of BL Lac objects. 

Although a tight correlation between Lrsq and is found for 7-ray BL Lac objects, 
we flnd that only a small fraction of BL Lac objects were detected by EGRET, which may 
imply that the detection of blazars in 7-ray energy band is not only dependent on their radio 
emission, but also on their 7-ray activity. The 7-ray emission is closely related with radio 
emission for BL lac objects based on SSC models. More distant or/and fainter BL Lac objects 
can be detected if the sensitivity is improved. Abdo et al. (2009a) found that the redshift 
distribution of Fermi BL Lac objects {z < 1.2) is similar to that of the DXRBS BL Lac sample 
{z ^ 1.0) (Padovani et al. 2007). This may imply that the radio-gamma correlation found in 
this work has not been affected by the redshift distribution of 7-ray BL Lac objects. Thus, 
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we induce a 7-ray duty cycle 5^ for BL Lac objects, which can be estimated from the ratio of 
7-ray BL Lac objects to the total (e.g., Stecker & Salamon 1996; Vercellone et al. 2004). There 
are 36 BL Lac objects were detected by EGRET and Fermi, among which the lowest radio flux 
density at 5 GHz is ~0.25 Jy. The number of BL Lac objects above this flux density limit with 
the same sky coverage of the survey carried out with EGRET/Fermi is ~330 (e.g., see figure 4 
in Padovani et al. 2007), which leads to 5^ ~ 0.11. Vercellone et al. (2004) estimate the 7-ray 
activity using three different methods. Our result is consistent with the maximum distribution 
of blazars active fraction computed in Vercellone et al. (2004) (see Fig. 14 of Vercellone et al. 
(2004)). 

The 7-ray LF of BL Lac objects can then be calculated with 

$,(L„.) = 5,<|.r(Lr,5G,.)^^ = ^^|^$r(V5G,^), (2) 

where $^ is the 7-ray LF, and the radio LF $p{ is given by Padovani et al. (2007). We can use 
a simple GLF model defined as 

$^(L^,2) cxL^, (3) 

In figure 2, we compared our 7-ray LF (/3 = —2.13) derived from the radio LF and that directly 
derived with a sample of 7-ray BL Lac objects detected by Fermi (Abdo et al. 2009a). The 
slope of our 7-ray LF is well in agreement with the value of —2.17 ±0.05 reported for Fermi 
BL Lac objects. Our 7-ray LF is slightly less than that of (Abdo et al. 2009a), but taking into 
account the observational error, we can also say our 7-ray LF is roughly consistent with that 
of Fermi. 

The number count of BL Lac objects as a function of 7-ray flux above 100 MeV can be 
calculated with the derived 7-ray LF, 

/■Zm fly f 

iV(>0=/ -rd^ / ULua.z)dL,,,. (4) 



dz 

where = 1 is adopted in our calculations, because the radio LF for BL Lac objects was 
derived with a sample of sources with redshifts ^ = — 1 (see Padovani et al. 2007, for the 
details). This means that the results derived here are only the lower limits. In fact, almost all 
7-ray BL Lac objects detected have redshifts (only 2 of 40 have redshifts slightly higher 

than unity). We plot the number count of 7-ray BL Lac objects derived from the LF calculated 
with equation (2) in figure 3. Obviously, the increasing number of sources can be detected 
following a decreasing flux limit, such as ~ 500 at /umit = 0.06/i^f^^ (the Fermi sensitivity 
for one year) (Gehrels & Michelson 1999), and A^ ~ 1000 at /umit = 1/30/if^f-^^ (the Fermi 
sensitivity for two year) (Atwood et al. 2009). 

The EGRB contributed by all BL Lac objects can be calculated with the 7-ray LF 
derived with equation (2), 
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where the average photon spectral index T = 2.04 for BL Lac objects is adopted, which is 
shghtly higher than F = 1.99 ±0.22 given by Abdo et al. (2009a). This is due to that our 
sample includes 9 EGRET sources, of which the photon spectral indexes are all greater than 
2.0 except 1011+496. The contribution of the BL Lac objects to the EGRB as a function of 
7-ray flux limit /umit is plotted in figure 4. 

4. Discussion 

It was suggested that the 7-ray radiative mechanisms are different for quasars and BL 
Lac objects, i.e., the EC mechanism may be responsible for quasars, while the SSC is for BL Lac 
objects (e.g., Dondi & Ghisellini 1995). A linear relation between radio and 7-ray emission is 
expected for AGNs, if the SSC mechanism is responsible for 7-ray emission (see, e.g., equation 
(28) of Dermer et al. 1997). The correlation between radio and 7-ray emission for BL Lac 
objects found in this work is very close to a linear one as described in equation(l), which seems 
to support the SSC mechanism for 7-ray emission from BL Lac objects. 

The variability of 7-ray emission in 10 sources of 42 Fermi BL Lac objects have been 
detected. The faction is much less than that of FSRQs (45/57) (Abdo et al. 2009a), which 
implies that the variability of 7-ray emission from BL Lac objects could be weaker than that 
from FSRQs. The detection rate in 7-ray band may be dominantly related to the duty cycle 
of 7-ray BL Lac objects. Vercellone et al. (2004) investigated the 7-ray activity of EGRET 
blazars and estimated their 7-ray duty cycle. They found that about 48 percent of the sources 
fall into 6^ = — 5%, and 74 percent sources are in 5^ < 10%. The duty cycle 5^ ~ 0.11 derived 
in this work is roughly consistent with the results in Vercellone et al. (2004) (while their 7-ray 
duty-cycle was computed without a clear distinction between BL Lacs and FSRQs.). 

In this work, the 7-ray LF of BL Lac objects is derived from the radio LF by using the 
empirical correlation of 7-ray emission with radio emission, which can almost reproduce the 
7-ray LF calculated directly with the 7-ray BL Lac objects detected by Fermi/EGRET (see 
figure 2). Based on the derived 7-ray LF for BL Lac objects, we can calculate the number 
count of 7-ray BL Lac objects as a function of flux limit in 7-ray band. We find that about 
1000 BL Lac objects will be detected by Fermi, if its sensitivity is 30 times higher than that of 
EGRET at 100 MeV (Gehrels & Michelson 1999). Cao & Bai (2008) 's estimate show that about 
1200 quasars can be detected by Fermi based on the EC mechanism for their 7-ray emission. 
The recently released bright AGN sample detected by Fermi together with the sources already 
detected by EGRET with high- confidence leads to a sample of 135 7-ray blazars (77 quasars, 
51 BL Lac objects and 7 blazars with uncertain classification) . The ratio of quasars to BL 
Lac objects for the present sample (LBAS+EGRET) is about 1.5, which is slightly higher than 
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Fig. 3. The number count of 7-ray BL Lac objects as a function of flux limit at 100 MeV. The dotted 
line and dot-dashed line represents the sensitivity of Fermi. We note /min ~ 5 x 10~^ ph cm~^ for 
EGRET (see, e.g., Cao & Bai 2008). The flux limit of Fermi at 100 MeV can be 30 times lower than that 
of EGRET for a two year all-sky survey (Gehrels & Michelson 1999) and the sensitivity for one year will 
be 3.0 X W~^ph cm~'^ s~^{Atwood et al. 2009). 
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Fig. 4. The contribution of the BL Lac objects with 7-ray flux / > /nmit- The dashed line indicates the 
measured EGRB with photon energy above 100 MeV. We use the value of the EGRB obtained by Strong 
et al. (2004), (1.14 ±0.12) x 10~^photons cm~^ s~^ sr~^. The symbols of dotted line and dot-dashed line 
as in Figure 3. 



1200/1000, which imphes that more BL Lac objects are expected to be detected by the Fermi 
in the future. 

Integrating the derived 7-ray LF, we calculate the contribution of all BL Lac objects 
to the EGRB, which accounts for ~45% of the EGRB (see figure 4). Cao & Bai (2008) 's 
calculation show that all radio quasars and FR II radio galaxies contribute about ~32% of 
the EGRB. This means that the contribution of BL Lac objects is similar to that of quasars, 
which is roughly consistent with the predicted number ratio of quasars to BL Lac objects to 
be detected by Fermi. The contribution of all blazars (quasars+BL Lac objects) can account 
for ~77% of the EGRB, which implies there still some space (^ 23%) is left for other sources 
to the EGRB. Considering that our calculations are limited to z = ~ 1 for BL Lac objects, 
the contribution of the sources other than blazars to the EGRB should be lower than the value 
derived in this work. We also noticed the detection of 7-ray emission from relativistic jets in 
the narrow-hne Seyfert-1 galaxies (PMN 0948+0022) (Abdo et al. 2009b), which means that the 
mechanism of 7-ray radiation from narrow-line Seyfert-1 galaxies is similar to that for blazars. 
Considering that a fraction of narrow-line Seyfert-1 galaxies are radio-loud (e.g. Komossa et 
al. 2006; Yuan et al. 2008) and some of them show evidence for relativistic jets (e.g. Gu & 
Chen 2010), the contribution of radio-loud narrow-line Seyfert-1 galaxies to the EGRB may be 
important, however, the detailed calculation is beyond the scope of this work. 

We thank the referee for his/her helpful comments, and J. M. Bai for helpful discussion. 
This work is supported by the NSFC (grants 10773020, 10821302 and 10833002), the CAS 
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Table 1. Literature Data for 7-ray BL Lac objects 
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